Analysis of nonperturbative fluctuations in a triple-well potential 
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We consider the quantum tunneling phenomenon in a well-behaved triple- 
well potential. As required by the semiclassical approximation we take into 
account the quadratic fluctuations over the instanton which represents as 
usual the localised finite-action solution of the euclidean equation of motion. 
The determinants of the quadratic differential operators at issue are evaluated 
by means of the Gelfang-Yaglom method. In doing so the explicit computation 
of the conventional ratio of determinants takes as reference the harmonic 
oscillator whose frequency is the average of the individual frequencies derived 
from the non-equivalent minima of the potential. Eventually the physical 
effects of the multi-instanton configurations are included in this approach. As 
a matter of fact we obtain information about the energies of the ground-state 
and the two first excited levels of the discrete spectrum at issue. 
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The semiclassical analysis of the tunnehng phenomenon, ranging from periodic-potentials 
in quantum mechanics to Yang-Mills models in field theories, has been performed in a 
systematic way by means of the instanton method. To be precise, the instanton represents 
itself a localised finite-action solution of the euclidean equation of motion. In principle 
one finds the classical configuration and subsequently evaluate the quadratic fluctuations. 
The functional integration at issue is solved in terms of the gaussian scheme except for the 
zero-modes which appear by virtue of the translational invariance of the model. Once the 
collective coordinates technique is incorporated the gaussian integration takes into account 
the directions orthogonal to the zero-modes. As a functional determinant includes an infinite 
product of eigenvalues, the result should be a highly divergent expression. However we can 
regularize the fluctuation factors by means of the ratio of determinants. 

Let us describe in brief the instanton calculus for the one- dimensional particle as can 
be found for instance in [|l|. Unless otherwise noted we assume that our particle moves 
under the action of a confining potential V{x) which yields a pure discrete spectrum of 
energy eigenvalues. In addition we choose the origin of the energy so that the minima of 
the potential satisfy V{x) = 0. If the particle is located at the initial time ti = — T/2 at 
the point Xj while one finds it when tf = T/2 at the point Xj, the functional version of the 
non-relativistic quantum mechanics allows us to express the transition amplitude in terms 
of a sum over all paths joining the world points with coordinates (—T/2, Xi) and (T/2, Xf). 
Making the change t —iT, known in the literature as the Wick rotation, the euclidean 
formulation of the path-integral reads 



where H represents the hamiltonian of the model, the factor N{T) serves to normalize the 
amplitude conveniently while [dx] indicates the integration over all functions which fulfil the 
boundary conditions at issue. In addition we resort to the euclidean action Se, i.e. 




(1) 




dr 



(2) 



2 



whenever the mass of the particle is set equal to unity for notational simplicity. This 
semiclassical approximation takes for granted the existence of a Xc{t) configuration which 
represents a stationary point of the euclidean action. Next we perform the expansion of a 
generic x{t) according to 



j 

where Xj (t) stand for a complete set of orthonormal functions 

/ Xo(t) Xkir) dr = Sj 

J-T/2 



(3) 



(4) 



vanishing at our boundary, i.e. Xj{i:T/2) — 0. As a matter of fact the eigenfunctions Xj(t) 
(with eigenvalues ej) appear associated with the so-called stability equation given by 

d^Xjir) 



+ V"[x,ir)]xjir) = e,x,{r) 



(5) 



In doing so the euclidean transition amplitude reduces itself to 



< Xf\eM-HT)\xi >= N{T) eMSeo) 



-1/2 



(6) 



where Seo represents the classical action associated with the configuration Xc{t) while the 
product of eigenvalues is usually written as 



Det 



dr^ 



+ V"[x.{t)] 



-1/2 



(7) 



according to a notation which obviously has its origin in the finite-dimensional case. To fix 
the factor N{T) we resort to a well-known problem where 



V{x) ^-x' 



(8) 



so that V"{x — 0) — u'^. As corresponds to the harmonic oscillator the relevant amphtude 
is given by 



<Xf^O\ ex.-p{-HhoT)\xi = >= N(T) i^Det 



dr^ 



-1/2 



(9) 



Next the evaluation of @ is possible according to the general method exposed in 0]. To 
sum up 

/jy\l/2 

< X/ = 0| exp{-HhoT)\xi = >= i-j (2 sinh z/T)"^/^ (10) 

In the following we consider the triple- well potential V{x) given by 

V{x) = ^x^ix^ - If (11) 
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From a classical point of view we find three minima located at X- = —1, x+ = 1 and 
Xo = 0. Two of them, namely x_ and equivalent since can be connected by means of 

the discrete symmetry x — x at issue. However the third minimum Xo is invariant under 
the action of such a symmetry. In other words, the central vacuum is not identical with the 
others located at both sides. When considering the limit uj"^ ^ 1 the energy barriers are high 
enough to decompose the system into a sum of independent harmonic oscillators. However 
the existence of finite barriers between the different wells of the potential yields a relevant 
tunneling phenomenon so that ultimately the symmetry x — > —x is not spontaneously broken 
at quantum level and the expectation value of the coordinate x computed for the ground- 
state is zero as corresponds to the even character of the potential V{x). The question we 
wish to address now is the explicit description of the tunneling in the euclidean version of the 
path-integral. As regards the one-instanton amplitude we evaluate the transition amplitude 
between Xo = and x+ = 1. In doing so we need a classical configuration with Xj = at 
ti = —T/2 while Xf = 1 when tf = T/2. It is customarily assumed that T — oo mainly 
because the explicit solution of the problem is much more complicated for finite T (more on 
this point later). 

As a matter of fact we can find the explicit form of the instanton Xc{t) just by integration 
of a first-order differential equation, i.e. 

^ = ±,/m^) (12) 

where we recognize the quantum mechanical version of the Bogomol'nyi condition |^. Now 
we solve (0) by a simple quadrature so that 
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Xr\r 



[1 + tanhu;(r - Tc)]/2 



(13) 



where the parameter Tc indicates the point at which the instanton makes the jump. As 
expected equivalent solutions are obtained by means of the transformation r — r and 
2^c(t) — ^ —Xc{t) so that adjoint minima of the potential can be connected by means of a 
topological solution. In addition we have that Seo = 00/ 4. On the other hand we need 
classical configurations for which x_ = and a;+ = 1 at large but finite values r = ±T/2. 
However the explicit form of the instantons that appear in the literature corresponds to 
infinite euclidean time. Fortunately the difference is so small that can be ignored mainly 
because we are interested in the limit T — ^ 00. In any case our description of the one- 
instanton amplitude between x_ = and = 1 takes over 



<Xf = l\ exp(-iJT)|x, = >= N{T) j^et 



-1/2 



Det [~{dydT^) + V"[xc{T)]] 



-1/2 



exp(-5'eo) 



(14) 



t Det[-{dydT^) + iy^] J 
where we have multiplied and divided by the determinant associated with the harmonic 
oscillator of frequency u. As regards the determinant built over the instanton itself we find 
a zero-mode Xo{t) which could jeopardize the computation process. However this eigenvalue 
€0 = comes by no surprise since it reflects the translational invariance of the system. As a 
matter of fact one can discover the existence of a zero-mode starting from (0). Including 
the adequate normalization one can check that 

1 dXr 



Xo{r) 



(15) 



/ Seo dT 

is just the solution of (|^) with eo = 0. The way out of this apparent cul-de-sac is simple. 
The integration over Cq (see (P)) becomes equivalent to the integration over the center of the 
instanton t^. To fix the jacobian of the transformation involved we take a first change such 
that 



Ax(r) = Xo{t) Aco 



(16) 



According to the general expression written in we find that under a shift Atc the 
effect should be 

Ax(r) = -/S^ Xoir)Ar, (17) 

Now the identification between (^) and (|T^ yields 

dco = \/sZ (18) 

where the minus sign disappears since what matters is the modulus of the jacobian at issue. 
In doing so we have that 

[ Pet [-(rfVrfr^) + V"[x,{t)]] \ _ 
j Det[-{dydT^) + u^] J 

i Det'[-{dydr') + V"[x^{r)]] \-'/' 

\ Det [-{dydr^) + / V 27r ^ 

where Det' stands for the so-called reduced determinant once the zero-mode has been re- 
moved. In order to make an explicit computation of the quotient of determinants we resort 
to the well-grounded Gelfand-Yaglom method so that only the knowledge of the large-r 
behaviour of the classical solution Xc(r) is necessary [If. Being O and P a couple of se- 
cond order differential operators, whose eigenfunctions vanish at the boundary, the quotient 
of determinants is given in terms of the respective zero-energy solutions /o(t) and 5'o(t) 
according to 

DetO fo{T/2) 



DetP goiT/2) 
whenever the eigenfunctions at issue satisfy the initial conditions 



(20) 



/„(-T/2) = (?,(-T/2) = 0, ^(-T/2) = ^(-T/2) = 1 (21) 

As regards the zero- mode 5'o(t) associated with the harmonic oscillator of frequency u 
we find 

^7o(r) = ^ sinh[z/(r + T/2)] (22) 



so that now we need the exphcit form of the solution /o(t) which corresponds to the topo- 
logical configuration written in (p^. Starting from the Xo{t) zero- mode we can construct a 
second solution 2/0 (t) given by 

yo{r)=x^{r) (23) 
Jo x^[s) 

Accordingly we may summarize the asymptotic behaviour of Xo(r) and VoiT) as follows 

{Cexp(— 2to'r) if r ^ 00 
(24) 
-Dexp(ci;r) if r — > — 00 

{exp{2LJT) / AujC if r ^ 00 
(25) 
— exp^—cor) /2ujD if r ^ —00 

where the constants C and D derive from the explicit form of the derivative of (|T3|). We 

proceed to investigate the particular solution /o(t) which is the one we are really interested 

in. Starting from the linear combination of Xo{t) and ?/o(t) given by 

fo{T) = Axo{t) + By,{T) (26) 

the incorporation of the initial conditions at issue leads us to 

fair) = x,(-T/2)y,(r) - y„(-T/2)a;,(r) (27) 

From this expression, which is exact, we can extract the asymptotic behaviour of /o(t), 

i.e. 

/o(T/2) ~ ^ exp{u;T/2) tf T ^ 00 (28) 

Next we need to take into account the lowest eigenvalue of the stability equation to obtain 
the right value of the quotient of determinants. From a physical point of view we can explain 
the situation as follows: the derivative of the topological solution does not quite satisfy the 
boundary conditions for the interval (— T/2,T/2). When enforcing such a behaviour, the 
eigenstate is compressed and the energy shifted slightly upwards. In doing so the zero-mode 
Xo(t) is substituted for the /a(t) which corresponds to 



- + V"[x^{T)]h{r) = Xhir) (29) 

whenever 

/a(-T/2) = /a(T/2) = (30) 
To lowest order in perturbation theory we get 

h{r) ~ fair) + A ^ (31) 

so that ultimately 



A=0 



fx{r) = Mr) + A / [x,{T)y,{s) - yo{T)x,{s)] Us) ds (32) 

J-T/2 

The asymptotic behaviour of f°{T), Xo(r) and Voij), together with the condition 
f\(T/2) = allows us to find the lowest eigenvalue A, i.e. 

\ = 2uDexp{-ujT) (33) 

In such a case the Gelfand-Yaglom method provides us with the final expression of the 
quotient of determinants whenever we choose for the frequency u of the harmonic oscillator 
of reference the average of the frequencies of the central and lateral wells. In other words 
u = 3uj/2. Notice the difference with the well-grounded double-well model where the two 
minima of the potential are equivalent so that the aforementioned average is not necessary. 
Going back to the explicit form of Xc{t) (see (p!3|)) we find C = D = uj/y/S^. Armed with 
this information we can write the one-instanton amplitude between the points Xi = and 
Xf = 1 in terms of 

<Xf = l\exp{-HT)\xi = >= 
3^x1/2 FY' 

— (2sinh3a;T/2)-^/' J \ ^MSeo) uj dr, (34) 

Apart from the first factor, which represents the contribution of the harmonic oscillator 
of reference, we get a transition amplitude just depending on the point Tc at which the 
instanton makes precisely the jump. 
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Although all the above calculations were carried out over a single instanton, it remains 
to identify the contributions which take into account the effect of a string of widely sepa- 
rated instantons and antiinstantons along the r axis. It is customarily assumed that these 
combinations of topological solutions represent no strong deviations of the trajectories just 
derived from the euclidean equation of motion without any kind of approximation. We shall 
compute the functional integral by summing over all such configurations, with k instantons 
and antiinstantons centered at points Ti,...,Tk whenever 

T T 

- - < n < ... < Tfc < - (35) 

Being narrow enough the regions where the instantons (antiinstantons) make the jump, 
the action of the proposed path is almost extremal. We can carry things further and assume 
that the action of the string of instantons and antiinstantons is given by the sum of the k 
individual actions. This scheme is well-known in the literature where it appears with the 
name of dilute gas approximation In addition the translational degrees of freedom of 
the separated k instantons and antiinstantons yield an integral of the form 



rT/2 rTk rT2 (ujTY 

I udTj / LUdTk-l-.. / LudTi = — — — (36) 

J-T/2 J~T/2 J-T/2 k\ 



As regards the quadratic fluctuations around the k topological solutions we have now 
that the single ratio of determinants transforms into 

r^'^V' f2sinh3c.T/2)-/^ [ Pef [-{dVdr^) + V"[x.{r)]] \-"' 
(^-j (2smh3c.T/2) Det\-{dVdr^) + ~" 



— y^'exp(-3cur/4) 



/ Pet' [-{dydr^) + V"[xc{t)]] 
{ Det[-{dydT^) +9ujy A] 



-1/2] ^ 



(37) 



,27r. 

according to the limit of the factor associated with the harmonic oscillator when T is large. 
In addition, when going to the dilute-gas approximation we must consider a set of instantons 
and anti-instantons so that each topological configuration starts where its predecessor ends. 
For the amplitude at issue the total number k of instantons plus anti-instantons must be 
odd. On the other hand, the combinatorial factor F associated with the number of possible 



configurations corresponds to F = 2*^*^"^^/^ since the closed paths starting and coming back 
to the point Xo = require in a systematic way instanton-anti-instanton pairs. Notice the 
difference with the double-well potential where the instantons strictly alternate with the 
anti-instantons since the problem has only two minima so that the combinatorial factor is 
not necessary at all. Next we can write the complete transition amplitudes for the triple- well 
potential so that 

<xj = l\ eM-HT)\x. = >= (-) exp(-3.;T/4) E ^(^j + 1)! ^^^^ 

where d stands for the so-called instanton density, i.e. 

/ g . — 

\ nr~ y^eo exp{-Seo) (39) 

V OTT 

In summary 

/Q^N 1/2 

<Xf = l\exp{-HT)\xi = >= — exp(-3c<;T/4) smh{uTd) (40) 



Air, 

Going back to the formal euclidean transition amplitude written in(0), the insertion of 
the pure discrete spectrum of energy eigenfunctions, namely 

H\n >= En\n > (41) 

allows us to write that 

< Xf = 1\ exp{~HT)\xi = >= ^exp{-EnT) < Xf = l\n X n\xi = > (42) 

n 

Let us denote by Eg, Ei and E2 the energies of the ground-state and the two first excited 
levels of our problem. As the triple-well potential we are dealing with is even we know that 
< l|xj = > vanishes so that the limit T — >^ 00 in (|40|) provides us with the following energy 
eigenvalues 

SciJ , ^ 307 ^ Scu , 
Eo = — -ujd, ^1 = E2 = — + ujd (43) 

In other words, the tunneling phenomenon transfers along the whole real axis the gaussian 
wave functions constructed on the points x_ = —1, x+ = 1 and Xo = 0. As expected the 
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splitting term results proportional to the barrier-penetration factor. The average of the 
harmonic frequencies over the non-equivalent minima of the potential serves as the central 
position for the splitting. 
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